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Abstract. It was recently shown that the point spectrum of the separated Coulomb- 
Dirac operator Ho{k) is the limit of the point spectrum of the Dirac operator with 
anomalous magnetic moment Ha{k) as the anomaly parameter tends to 0; this spectral 
stability holds for all Coulomb coupling constants c for which Ho{k) has a distinguished 
self-adjoint extension if the angular momentum quantum number k is negative, but for 
positive k there are certain exceptional values for c. Here we obtain an explicit formula 
for these exceptional values. In particular, it implies spectral stability for the three- 
dimensional Coulomb-Dirac operator if |c| < 1, covering all physically relevant cases. 



1. Introduction 



By separation of variables in spherical polar coordinates, the Dirac operator of rela- 
tivistic quantum mechanics with a Coulomb potential, 

Uo := -ia ■ V -f /3 + 

where c < and ai, a2, eta and /3 = cto are symmetric 4x4 matrices satisfying the 
anticommutation relations 

aiUj + ajai = Sij (i, j G {0, 1, 2, 3}), 

is unitarily equivalent to a direct sum of one- dimensional Dirac operators on the half- 
line, 

d /t c 

Hoik) = -ia2-r + (T3 + -c^i + - (re (0, cx))), 
dr J. J. 

where o)'^^~(^i 0^) ^"^ ~ ^1 ) Pauli matrices 

and k & TL \ {0} is the angular momentum quantum number [13, Appendix to Section 
1]. For an electron orbiting a nucleus of charge number Z G N, the Coulomb coupling 
constant is c = —Za, with the Sommerfeld fine structure constant a ~ 1/137. 
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HQ{k) is essentially self-adjoint on its minimal domain if and only if < fc^ — 1/4 
[13, Theorem 6.9]. For G (/c^ — 1/4, /c^), there still exists a distinguished self-adjoint 
extension, characterised by the requirement that the wave-functions in its domain have 
a finite potential (or kinetic) energy [13, Theorem 6.10]. 

Pauli suggested a modification of the Dirac operator which takes into account the 
anomalous magnetic moment of the electron (for the historical background cf. [7]). 
With suitably normalised constants, the operator with a potential V for an electron of 
magnetic moment (1 -|- a)iJ,B, where //s is the Bohr magneton, takes the form 

Ha = -ia ■ V + P + V - iaa ■ VF. 

In the case of the Coulomb potential, the corresponding half-line operators after variable 
separation will be 

Ha{k) = -i(J2:^ + (73 + f - + 4 Vl + - (^ ^ (0' '^))- 

dr \r J r 

The mathematical investigation of the properties of this operator was initiated by 
Behncke [2, 3, 4]. Ha{k) is essentially self-adjoint on its minimal domain for all values 
c < of the coupling constant (cf. [1, 5]). The essential spectrum of Ha{k), as that of 
Ho{k), is (-00, -1] U [1, oo). 

Let us now focus on the case < k^, so that a distinguished self- adjoint realisation 
of Ho{k) exists. We always assume a < in the following; the case of positive a can be 
reduced to this (with a sign change of A;) by a suitable unitary transformation. Although 
only integer values of k are relevant for the three-dimensional operator, we admit general 
non-zero real values of k. As is well known, -ffo(^) has infinitely many eigenvalues in 
the spectral gap (—1, 1), which accumulate at 1. One would expect that the eigenvalues 
of Ha{k) will be perturbations of those of Ho{k), such that each eigenvalue of HQ{k) 
will be the limit of exactly one eigenvalue branch of Haik) in the limit a — > {spectral 
stability, cf. [8, Chapter VIII §1.4]). This expectation is partly corroborated by the 
strong resolvent convergence of Ha{k) to Hoik), at least for < P - 1/4, which 
implies that the spectrum cannot suddenly expand in the limit. Nevertheless, due to 
the strong singularity of the a/r^ term at the origin, this limit problem is highly non- 
trivial and does not subject itself to the standard perturbation techniques. Behncke [4, 
Theorem 2] was able to prove spectral stability for < A;^ — (3/2)^ if A; < 0, and for 
c2 < A; - 5/2 if A; > 0. 

The surprising asymmetry with respect to the sign of k is not an artefact of 
Behncke's method, but inherent in the problem, as observed in the recent treatment [7], 
where spectral stability was proved, based on an asymptotic analysis of the (non-linear) 
Priifer and Riccati equations equivalent to the Dirac system 

{Ha{k) - \)u = 0, (1) 

for aU c e (A;, 0) if A; < 0, and for aU c e (-A;, 0) \ {cq, ci, . . .} if A; > 0. Here cj e {-k, 0), 
with Cj > Cj+i, are certain exceptional values of the coupling constant at which a 
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transition in the behaviour of the eigenfunctions of Ha{k) occurs: for c e (c^, c„_i), 
the eigenfunctions of Ha{k) show n additional rapid oscillations very close to the origin, 
compared to the corresponding eigenfunctions of Ho{k). It remains a fairly subtle open 
question whether or not spectral stability holds if c is equal to one of the exceptional 
values. 

The present note is devoted to a study of the number and positions of the excep- 
tional values of c for a given k > 0. 



2. Exceptional values 



The exceptional values for c appear in a stability analysis of the differential equation 
[7, Eq.(l)] 

g'd'ig) = c + (/c - -) sin 2'd {g > 0). (2) 
Q 

This equation arises from the Priifer transformation (cf. Appendix A) of the Dirac 
system (1) after omitting the lowest-order term at 0, us — A, and rescaling r = \a\g to 
absorb the parameter o. 

In the limit g — > cxo, the right-hand side of (2) — which is 7r-periodic in 7? — has 
asymptotic zeros 'i?±(c, k) satisfying < ^?_(c, /c) < 7r/4 < 'i?+(c, k) < 7r/2, 



■ r^n/7\ c n/7\ Vk^ — 

sm 2v±(c, k) = — — , tan v±{c, k) = . 

k —c 

An asymptotic study of the direction field of (2) for g and ^ — > oo yields the 
following result, which shows that this differential equation has (up to addition of a 
constant multiple of tt) a unique unstable solution at and at oo (see [7, Lemma 2.1, 
2.2]). 



Lemma 1. Let k > and c e (— /c, 0). Then the following statements hold. 

a) There is a unique solution 'do{ ■ , c) of (2) such that lim 'do{g, c) = tt. 

All other solutions either differ from -do by a constant integer multiple ofir, or else 
satisfy lirn ■i?(^, c) = ^ mod tt. 

For fixed ^ > 0, '&o{g^ ■ ) is continuous non- decreasing. 

b) There is a unique solution '&oo{ ■ , c) of (2) such that lim 'dooiQ, c) = ''?-(c, k). 

Q—^OO 

All other solutions either differ from t^oo by o constant integer multiple of tt, or 
else satisfy lim z?(^, c) = 'i?+(c, k) mod n. 

Q — >00 

For fixed g > ''9oo{q, ■ ) is continuous and strictly decreasing. 



The solution ■do{ ■ , c) asymptotically corresponds to the Priifer angle of the L^(0, oo) 
solution of (1). 

The exceptional values are defined as those values of c for which 'do{ • , c) and 
■&oo{ ■ 1 c) match up mod tt, so that the unstable solution of (2) at is unstable at 
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infinity as well. More precisely, in view of the monotonicity properties of t^q and "^^q 
with respect to c, we have 

lim 'do{g, Cm) = ^^-(cm, k) - rmr (m = 0, 1, 2, . . .). 

By the transformation g = e*, = <^(log g), the equation (2) is equivalent to 

ip'{t) = c + {k - e"*) sin 2ip{t) {t e R), (3) 

the differential equation for the Priifer angle (p (cf. Appendix A) of a R^-valued solution 
u of the Dirac system 

+ (/e - e-')ai + c) u{t) = 0, (4) 

which can be rewritten using the definition of the Pauli matrices as 

u'l = (e~* — k)ui — CU2, 
u'2 = cui -\- {k — e~^)u2- 

In (4), c takes on the role of a spectral parameter, while the coefficient of ai can be 
interpreted as a constant mass term with an exponentially decaying perturbation. In 
the sense of the analogue of Kneser's Theorem for Dirac systems [9], this perturbation 
is subcritical, and indeed the method developed in [9] can be used to show that, for each 
/c > 0, there are only finitely many exceptional values for c. Moreover, an asymptotic 
analysis of (4) along the lines of [10; see also 12] reveals that the number of exceptional 
values is asymptotic to k in the limit A; — > 00 with fixed ratio c/k. 

This was verified in a computational investigation of (4) based on a piecewise- 
constant approximation of the exponential function, following the approach of [11]. 
The numerical findings suggested that the number of exceptional values in (— /c, 0) is 
always equal to the unique integer in [fc — 1, /c); more precisely, whenever k reaches an 
integer value, an additional exceptional value, initially equal to — /c, appears and moves 
into the interval {—k, 0) as k increases. 

This regularity raised the suspicion that (4) has an underlying solvable structure, 
and indeed, on closer scrutiny it turns out that this equation can be analysed by means 
of a variant of the factorisation method, resulting not only in a proof of the above 
conjecture, but even in an explicit formula for the exceptional values, which makes any 
asymptotic and computational analysis of (4) superfluous. 

Theorem 1. For given A; > 0, the exceptional values for Ha{k) are given by 

Cn-i{k) = -\/2kn - (n e N, n < k). 

Consequently there are exactly N exceptional values in {—k, 0) if k & (N, N + 1], 

N e Kq. 
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In view of Theorems 1.1 and 1.2 of [7], this imphes in particular that the eigenvalues of 
Ha{k) converge to those of -ffo(^) as a ^ if < 2A; — 1, confirming Bchncke's linear 
bound (although with a different constant). However, his conjecture [4, p. 2558] that 
spectral stability may hold for all < /c^ — 5/2 turns out to be very questionable. 

Since the three-dimensional Coulomb-Dirac Hamiltonian Ha is the direct sum of 
Ha{k) for k e 7Z \ {0}, we can thus infer spectral stability for this operator in all cases 
in which Ho has a distinguished self-adjoint realisation. 



Corollary 1. Let Hq be the Dirac operator and Ha the Dirac operator with anoma- 
lous magnetic moment with a Coulomb potential with coupling constant c G (—1, 0). 
Then the eigenvalues ofHo are the limits of the eigenvalues of Ha as a — > 0. 



3. Proof of Theorem 1 



The following proof of Theorem 1 will be based on squaring the Dirac-type operator in 
(4), which is equivalent to deriving second-order differential equations for ui and U2. In 
a general situation, this is usually not a good idea except to obtain a quick and rough 
heuristic result, since it involves derivatives of the coefficients and consequently requires 
unnecessarily strong regularity. In the present instance, however, we are dealing with 
a specific Dirac system with analytic coefficients, and as it turns out, the resulting 
second-order equation system decouples and can be solved by the factorisation method 
in the relevant cases (the second-order equation (5) is closely related to the radial 
Schrodinger equation with a Morse potential [6, Section 5.2]). Here the factorisation is 
incidentally provided by the original Dirac system; see the discussion in Appendix B. 
The key observation is contained in the following result. 

Theorem 2. Let n > 0. Let vi{t, k.) := e~'^ ^-^t^ define Vj recursively for 
j + 1 by 

Vj+i{t, k) := {n + j- e-')vj{t) - v'j{t) (t e R). 
Then Vj{ ■ , k) is a nontrivial solution of 

v" = (e-^* - {2k + 2j - l)e-* + K^) v, (5) 
and it has the asymptotic properties 

lim Vjit, K) = 0, lim = — k, and lim —, -e =1 (7 e N). 



The proof can be done by induction with respect to j. The fact that Vj is a solution of 
(5) is readily verified by differentiating twice. The asymptotic properties are obvious 
for 

^^(^' ^) = e-t _ ^. 
vi(t, k) 
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Now assume that j e N is such that the assertion is true. If Vj+i were trivial, 

this would imply that (up to multiplication by a constant) Vj(t, k) = e^'^'^^^^'^'^ oo 
{t oo), contradicting the first limit property of Vj. The asymptotic properties of Vj+i 
are easily checked using the identity 



K + j - e * - -f:{t, k) 



Corollary 2. Let k > 0, j G N and Vj be as in Theorem 1. Then 
hm = ^ lim = 0. 



This is a direct consequence of Theorem 1 in view of 

Vj{t, k) 1 



Furthermore, the following conclusion can be verified by a straightforward calculation. 
Corollary 3. Let Vj be defined as in Theorem 1, and let j & ¥i , k > j . Then 

is a nontrivial solution of the Dirac system (4) with c = —^Jlkj — p- The Priifer angle 
(f of u satisfies 

j 

lim taiiipCt) = — =, lim tan ip{t) = 0. 

Since, with c as in Corollary 3, 



, , , k - V/c^ - c2 k - Jk'^ - 2kj + p 
tan V- (c) = = — 



-c ^2kj - p V^kj -p' 

where we have used j < k, the Priifer angle </? in Corollary 3 corresponds via '&{g) = 
<^(log q) to a solution of (3) which coincides (up to addition of integer multiples of tt) 
with both '!?o( ■ ) c) and "dooi ■ ) c); hence such values of c are exceptional values. 
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We now conclude the proof of Theorem 1 by showing, based on the continuity 

and monotonicity properties of ^nd i}^, that all exceptional values are obtained in 
this way. Corollary 3 specifies the asymptotic behaviour of the relevant angle functions 
mod tt; a study of the zeros of Vj reveals their global behaviour as shown in Lemma 2 
below. Theorem 1 then follows in view of 

^{q, -V^kj - f) = ipj{\og g,k- j) {g>0,k>0,j e N, j < k), 

where (with arctan = 0) 

(Pj{ ■ , n) '■= arctan( — ' — ^\/2kj + p) + n. 
Vj+l{ ■ , k) 



Lemma 2. Let ac > 0. Then Vj{ ■ , k) has exactly j — 1 zeros. The associated angle 
function (pj satisfies 

lim 'Pjit, k) = TT 

t—^ — 00 

and 

lim ^j{t, k) = i?_(-V2acj + j2, K + j)- {j - l)7r {j e N). 

Proof. A look at the direction field of (3) with k = k + j shows that (pj{ ■ , k), being 
the unstable solution at both ±oo, is monotone decreasing. vi{ • , n) has no zeros, so 
(fi{-,K) has no zeros mod tt, and the assertion follows for j = I. 

To conclude the proof by induction, assume now that j G IN is such that the as- 
sertion is true. Then (fij{ ■ , k) — 7r/2 has exactly j zeros mod tt, so • , n) has 
exactly j zeros. Since the ■ , «;) is a non-trivial solution of a second-order equa- 
tion, Vj_^_i{ • , k) does not vanish at zeros of • , k). Consequently, Vj+i{ ■ , k) and 
Vj+2{ ■ , 1^) have no common zeros. Hence i^j+ii " 7 ^) has exactly j zeros mod tt, and 
therefore must converge to ■d-{—-\/2Kj + k + j) — {j — l)7r as t — > 00. 



Appendix 

A. The Priifer Transformation 
Consider a general Dirac system 

d 

(-ia2- h m(x)(T3 -I- l{x)ai + q{x))u{x) = (x e /) (6) 

ax 

on an interval / C R with locally integrable, real- valued coefficients m, / and q (we have 
absorbed the spectral parameter in the latter). If tt is a solution of (6), then so is its 
(component-wise) complex conjugate, so it is sufficient to study R^-valued solutions. 
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Since (6) is linear, a non-trivial solution will never take the value ^ J , so u{x) traces 

out an absolutely continuous curve in the punctured plane as x varies. 
Introducing polar coordinates in this plane by writing 

f \ I / M cos \ 

^ ' ' ^ ^ ' y sm -dix) J 

where the absolutely continuous function •& (the Priifer angle of u) is determined 
uniquely up to addition of an integer multiple of 27r — or indeed tt, since —u is again 
a solution of (6) — , a straightforward calculation yields the differential equation for -d 
{Priifer equation) 

i}' — m cos 2-?? + I sin 2-?? + q. 

This non-linear equation is equivalent to (6), because one can recover u from by 
noticing that (choosing xq G /) 

r 

\u{x)\ — \u{xq) \ exp / (m sin 2^} — I cos 2'd) {x E I). 



B. Dirac systems and the factorisation method 

Known to students of quantum mechanics as a trick to treat the Schrodinger equation 
for the harmonic oscillator, the factorisation method is in fact a tool of much wider 
scope for studying second-order differential equations (see the extensive discussion in 
[6]). It is based on the close link between the eigenvalue equations 

A*Au = Xu and AA*v = Xv, 

where A is typically a first-order differential operator: a solution u of the first equation 
gives rise to a solution v = Au of the second. 

Here we remark on how this method can serve to solve eigenvalue equations for 
certain one-dimensional Dirac operators, which in turn represent a factorisation of the 
Sturm-Liouville operators arising as their formal squares. 

Starting from a Dirac system (6), or equivalently 

u'l = —lui + {m — q)u2, 
U2 = {m + q)ui + lu2, 

with absolutely continuous coefficients /, m, q, we find that ui and U2 satisfy the 
second-order equations 

u'l = (m^ - q^ + - l')ui + {m! - q')u2, 
u'^ = {rr? - 9^ + + l')u2 + (m' + ^O^i- 
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This ordinary differential equation system decouples if m and q are constant. If we 
assume this in the following and set c := m — q, d := m + q and a := cd, we are 
dealing with the Dirac system 

u[ = —lUi + CU2 

u'2 = dui + lU2 
and corresponding second-order equations 

u'l = {a + l^ - (7) 
u'i = (a + ^2 + i')u^_ (8) 

These two Sturm-Liouville equations differ only in the sign oil' . Hence we can construct 
a chain of interlocking equation systems, along with special solutions, in the following 
way. 

Theorem 3. Let [a^jnen 0, sequence of real numbers with ai — 0, and let 
(^n)nG]N be a sequence of real-valued absolutely continuous functions on I such that 

ll+l + I'n+i + On+l =ll- ln + an (n G N). 

Let vi{x) := exp li {x e /), and define recursively for n e N + 1 

Then Vn is a solution of (7) with a = an, I = In, o-'f^d of (8) with a = an+i, I = In+i- 
Furthermore, for n e N + 1 and arbitrary £ R \ {0}, 

u := ( dr,^^ \ 

is a solution of the Dirac system 

u[ = -InUl + {an/dn)u2 (9) 
u'2 = dnUi + lnU2- 

Proof. Because of ai = 0, the first equation of the Dirac system (9) with n = 1 
decouples from the second and can easily be solved to obtain the formula for vi . As the 
first component of a solution of this Dirac system, Vi satisfies equation (7) with / = /i, 
a = ai. In view of the assumptions on the coefficients this equation is the same as 
(8) with I = 2, a = a2, the Sturm-Liouville equation satisfied by the second component 
of a solution of the Dirac system for n = 2. The corresponding first component of this 
solution can easily be calculated from the second equation of the Dirac system. It will 
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depend on 6,2^ but note that this dependence can be ehminated by normahsation as far 
as obtaining a solution V2 of the Sturm-Liouville equation (7) with I = I21 a = a2 Ss 
concerned. 

This process is then iterated to find the solutions for n G {3, 4, . . .}. 
Examples. 

1. The classical application of the factorisation method to the one-dimensional Schrodinger 
equation for a quantum-mechanical harmonic oscillator, 

—v"{x) + x'^v{x) — Af (x), 

is captured in the above scheme if we choose ln{x) := a; (n G N). The hypothesis 
of Theorem 3 is then satisfied with := — 2(n — 1) (n G N), and we obtain the 
well-known solutions 

Vn{x) = (-— + x)''-^e-^ 
with eigenvalues A„ := 2n — 1 (n G N). 

2. The case of equation (5) is treated by choosing ln{x) := k + {n — 1) — (n G W); 
we immediately find vi{x) — "'"'^^ and the recursion formula of Theorem 2. 
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